Scars of the Wigner function. 
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We propose a picture of Wigner function scars as a sequence of concentric rings along a two- 
dimensional surface inside a periodic orbit. This is verified for a two-dimensional plane that contains 
a classical orbit of a Hamiltonian system with two degrees of freedom. The orbit is hyperbolic and the 
classical Hamiltonian is "softly chaotic" at the energies considered. The stationary wave functions 
are the familiar mixture of scarred and random waves, but the spectral average of the Wigner 
functions in part of the plane is nearly that of a harmonic oscillator and individual states are also 
remarkably regular. These results are interpreted in terms of the semiclassical picture of chords 
and centres, which leads to a qualitative explanation of the interference effects that are manifest 
in the other region of the plane. The qualitative picture is robust with respect to a canonical 
transformation that bends the orbit plane. 
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the spectral Wigner function JlO| , p~T] | . indeed, the com- 
putations presented in this letter show that individual 
Wigner functions, as well as their energy average can os- 
cillate with large amplitudes deep inside the energy shell. 

At points x = (q, p) inside the energy shell, the semi- 
classical spectral Wigner function is jl0],[ll) : 



Sixteen years have passed since Heller |l[ detected scars 
of periodic orbits in individual eigenfunctions of chaotic 
systems. Explanations in terms of wave packets 0] or the 
semiclassical Green function |^,|| do predict an enhance- 
ment of intensity near the projection of a Bohr-quantizcd 
periodic orbit. However, such theories apply to a collec- 
tive superposition of states near this quantization con- 
dition, fn spite of some quite striking visual evidence, 
the issue of scarring for individual states, is confused by 
the fact that different branches of the same periodic or- 
bit may overlap so that their contributions interfere in 
the position space and this must be superposed on the 
expected random wave background for the chaotic state. 
This has led to the need to make quantitative assessments 
of scarring strength [|]||. We here show that the phase 
space picture can provide sharper qualitative evidence of 
the influence of a periodic orbit on its Wigner functions. 

The old conjecture of Berry and Voros ra] that the 
Wigner functions for eigenstates of chaotic Hamiltonians 
are concentrated on the corresponding classical energy 
shells is widely disseminated ffl, as it is compatible with 
Schnirelman's theorem ||. Nonetheless, no such restric- 
tion arises in the more recent semiclassical theory for the 
mixture of these states over a narrow energy window, i. e 



W(x,E,e) =^2Aj(x,E) e - £t ' /h cos 



(1) 

The sum is over all the trajectory segments on the i?-shell 
with endpoints, Xj±, centered on x, as sketched in Figjl]. 
The action is merely the symplectic area, Sj = <f p • dq, 
for the circuit taken along the orbit from Xj_ to x J+ and 
closed by the chord — £,-(x). The time of traversal for the 
stretch along the trajectory is tj and e is the width of the 
energy window over which we average individual Wigner 
functions W n (x), i.e., 



■W n (x) , (2) 



We 



where D is the dimension of the phase space (q, p) . 
shall not be concerned with the Maslov phase jj, nor 
with the amplitude Aj(x,E), except to note that these 
are purely classical quantities that vary smoothly with 
x inside the shell, as compared with the high frequency 
oscillations of the cosine factor. 
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The highly oscillatory nature of the spectral Wigner 
function inside the energy shell is mostly washed out on 
projection, i.e. for true probability densities. However, 
the oscillations are unavoidable in the study of interfer- 
ence effects. 

The energy shell itself is a caustic, because the chord 
£ — > as x approaches the shell. In this limit, the con- 
tributing trajectories either shrink to a point, or they 
are very close to a periodic orbit. The modification of 
(|l|) leads to the Berry theory of scars 0, futher refined 
in (!(]] , for evaluation points x close to the energy shell. 

The point here is that large open segments of periodic 
orbits can be constructed by adding multiple windings 
to a primitive, small segment of a periodic orbit. If con- 
ditions for phase coherence, to be discussed, are satis- 
fied, we can then obtain a scar deep inside the energy 
shell. This scar is located at the two-dimensional cen- 
tral surface constructed by the centers of all the chords 
with endpoints on the periodic orbit. For the case where 
the phase space dimension D = 2, the central surface 
is the phase space region enclosed by the convex hull of 
the orbit that here is the energy shell 1 . For D > 4, 
the two-dimensional central surfaces are still bounded by 
one-dimensional periodic orbits that lie within the higher 
dimensional energy shell. In the simplest case where the 
orbit is plane and convex, the central surface coincides 
with the part of the plane within the orbit just as in the 
case where D = 2. If the orbit is not plane, then the cen- 
tral surface will be more complicated and it may exhibit 
singularities and self-intersections. In general it will not 
coincide everywhere with the invariant surface formed by 
changing continuosly the energy of this periodic orbit. 

We emphasize that these scars of the Wigner function 
that correspond to a very recognizable oscillatory pattern 
have special localization properties in phase space of di- 
mension D > 4. Previous explorations of phase space 
representations of eigenstates have been too blunt to dis- 
cern these fine features. The work of Feingold et. al. 
p3[ involves a projection from a three-dimensional sec- 
tion, which is fine for the Husimi function, but washes 
away the details of an oscillatory Wigner function. On 
the other hand, Agam and Fishman |l4"| restrict their in- 
vestigation to the neighbourhood of the orbit and, hence, 
to the energy shell. This is just the edge of the central 
surface. 

In this letter we tested our prediction for the simplest 
case of an unstable periodic orbit lying in a two dimen- 
sional plane. Accidentally this plane is a symmetry plane, 
however, we have checked that the enhanced amplitude of 
the spectral Wigner function accompanies the distortion 
of the orbit due to a nonlinear canonical transformation. 



Since the Wigner function is not invariant in this case, 
we thus obtain an essentially new system, in which the 
new central surface is not a symmetry surface. Thus, our 
scarring effect cannot be attributed to a special property 
of the symmetry plane. This case is essentially different 
from recent studies of higher dimensional systems where 
the scars arise over marginally stable invariant planes 
that are indeed special ]l2| |. 

The chord structure for segments of a periodic orbit 
is sketched in Fig.||. This is similar to a system with 
D = 2, for which Berry showed that there is only 
one chord for most points x. However, we must now dis- 
tinguish between the chords £ in = — £ out and the sum in 
(|l|) includes each different winding of the periodic orbit 
for primitive orbit segments, "in" and "ouf, associated 
with the two respective chords. Evidently all these "in" 
and "out" contributions build up if the energy is close to 
one of the values for which this periodic orbit is Bohr- 
quantized. We show in jl6| that the condition for these 
two contributions to be in phase is the same as that for a 
maximal contribution to the Gutzwiller trace formula . 
Therefore, near the quantization energy, the sum of all 
the contributions in ([!]), owing to chords in the periodic 
orbit, can be approximated by an expression analogous 
to Berry's simplest semiclassical aproximation jl5],[l7| f° r 
the Wigner function in D = 2 systems, but now evalu- 
ated over points x in the two-dimensional central surface. 
Hence, in this case, we have an enhanced contribution to 
the spectral Wigner function with the phase Si n (x)/Ti 
and the Maslov phase 7 = 0. The successive contours 
Si n ( x ) ^constant thus determine rings of constant phase 
along the central surface. Therefore, we predict rings 
of positive and negative amplitude superimposed on the 
background of contributions from uncorrelated trajec- 
tory segments that also contribute to the spectral Wigner 
function. The edge of this system of rings along the orbit 
itself is the object of the Berry theory ||, though it does 
not deal with multiple windings. 

We have tested this prediction for the coupled nonlin- 
ear oscillator (D = 4) 

H(q,p) = ^^ + 0.05g?+^ 2 -|) 2 , (3) 

sometimes referred to as the "Nelson Hamiltonian" . The 
classical periodic orbits together with the topology of 
their families in the energy. vs. period plot has been stud- 
ied in |l8j and the wave functions and their scars were 
discussed in [|l9| . Evidently the plane q\ = p\ = is clas- 
sically invariant, so we obtain an isolated periodic orbit 
on this plane for each energy (the vertical family). This 



1 Note that if the orbit is convex the central surface is just 
the phase space region inside it. 
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is then, the simplest case where the central surface is flat 
- the central plane -. The restriction of the Hamiltonian 
to this plane defines a harmonic oscillator, so there is a 
single chord with endpoints in the periodic orbit (except 
for sign) for each point {p2, 92) inside the orbit, except at 
q-2=P2 = 0. 

We computed the wave functions (q|n) for the eigen- 
states of (||) within the range of energies 0.821 < E < 
0.836 taking % — 0.05 and using a basis of eigenfunctions 
4'm\qi)4'n2\q2 - qj/2) , in the same way as in Q. 
The Wigner function for each state was then calculated 
directly by the double integral over Q = (Qi, Q2): 

W «fa'P) = 7^)2 / d Q (q + Q/2|n)(n|q-Q/2)e- 4E ^. 

(4) 

In FigJH we display the average over Wigner functions 
for the energy window chosen, surrounding the 11 'th 
Bohr energy for the vertical orbit that has Maslov in- 
dex 3. This is more convenient than smoothing with the 
Lorentzian window in (||) and should produce essentially 
the same results. 

As an example, Fig.|] shows some of the individual 
Wigner functions in this window for x in the (92, P2) 
plane. Immediatly, we notice that all these Wigner 
functions are remarkably regular, roughly in the region 
Q2 < 0. Indeed, Fig.|] coincides in this region with the 
Wigner function of the restricted Hamiltonian, also with 
respect to the precise phase of the oscillations. The indi- 
vidual Wigner functions follow the same phase contours, 
but there is a varying phase shift. If we extrapolate the 
semiclassical theory to energy smoothings of the order of 
the energy spacing, the conclusion is that there are no 
orbit segments, other than the periodic orbit itself, con- 
tributing to the Wigner functions on this surface up to 
the Heisenberg time. This is more dramatic for the state 
n = 305 whose Wigner function on all the central surface 
is almost the same of that of the restricted Hamiltonian 
(see Fig.[|). It is important to note that only this state 
seems to have a visual scar of the vertical orbit in po- 
sition space, although superimposed to a random wave 
background. 

For qi > 0, we have found other orbit segments that 
do not lie on the central plane, but which do contribute 
to the spectral Wigner function on it. Indeed, this is 
the case for segments of both the symmetric periodic or- 
bits shown in Fig.|| projected onto position space. Need- 
less to say, the multiple windings of these periodic orbits 
could produce total contributions of the same order as 
the plane orbit, if they are nearly Bohr-quantized. This 
supplies a qualitative explanation for the existence of in- 
terference patterns to the left of the Wigner functions in 
Fig.^ and |]. Notice that the second of the orbits in Fig.^ 
reaches into a region with qi greater than is attained by 
the plane periodic orbit. So we account for Wigner func- 



tions reaching outside the energy shell (which coincides 
with the periodic orbit on the central plane). Note that 
these are only particular examples of the many symme- 
tric periodic orbits whose central surface intersects the 
central plane along lines with qi > [l8| , [l9| |. Therefore, 
it will be hard to make quantitative predictions in this re- 
gion in contrast to that where (72 < 0, for which no such 
orbits of short period were found. The enhanced scar 
pattern of the Wigner function of the state n = 305 over 
all the central plane evidently washes out other possible 
contributions. 

The central surface of the harmonic orbit considered in 
this work coincides with an invariant plane of the Hamil- 
tonian (H) which is a reflection symmetry plane of the 
system. This might induce to the conclusion that the 
symmetry is the actual responsible for the scarring ob- 
served. We show that this is not the case by applying 
a nonlinear canonical transformation, described in [fl6|| , 
deforming the central plane q\ = p\ = into a new 
invariant surface (FigJ^ a). Now, the new central sur- 
face of the distorted periodic orbit no longer coincides 
everywhere with the invariant surface (FigJ^ b) losing 
in this way any symmetry. We calculated the Wigner 
functions for this new system, over the central surface of 
the periodic orbit, for eigenstates within the same range 
of energies surrounding the 11 'th Bohr level. Both the 
average and the individual Wigner functions, projected 
over the (92,^2) plane, display the same features as those 
shown in Fig.|| and Fig.|| respectively. Reference |l6| also 
considers the thickness of scar surfaces. 

Our main point is that individual trajectory segments 
do play a role in the spectral Wigner function. Indeed, 
phase coherent contributions from a periodic orbit can 
add up to scars in the form of concentric rings deep inside 
the energy shell. Eventhough the Berry- Voros hypothesis 
ties in more intuitively with the concept of quantum 
ergodicity, it should be noted that the chord picture of 
Wigner functions does not contradict Shnirclman's the- 
orem and subsequent exact results ||], because the os- 
cillations make negligible contributions to integrals with 
smooth functions. 

It would seem that the structure of interfering chords, 
from which the spectral Wigner function is built up, 
could only generate a messy picture where scars would 
be even harder to recognize than in wave functions for 
which quantitative methods are needed. We have now 
shown that the converse can be true in the simplest case, 
i.e. that the relatively low dimension of the central sur- 
face can allow it to miss, over an appreciable region, the 
contribution of chords from other orbits. In this region 
we obtain a regular ring structure, not only for the energy 
average, but even for the pure Wigner functions shown 
in Fig.|J. Though we cannot yet predict why the relative 
weight of this particular periodic orbit should vary so 
markedly among the states close to the Bohr-quantized 
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energy, we are now able to recognize in Fig.|4| a remark- 
able example of an individual scar. 
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FIG. 1. In general there are many orbit segments with end- 
points Xj_ and Xj+ centered on a given point x. The circuit, 
closed by the chord — £• = Xj_ — Xj_|_, defines the phase of the 
semiclassical contribution to (hi). 




FIG. 2. If the endpoints lie on a periodic orbit, we can add an 
infinite number of windings to the two shortest primitive segments, 
"in" and "ouf , divided by the chords £ in = — £ out . All the points 
x, centers of chords in the periodic orbit, form a two dimensional 
central surface. If the orbit is plane (as in the graph) the central 
surface is the part of the plane enclosed by the convex hull of the 
orbit. 
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FIG. 3. Averaged Wigner function on the (52, P2) plane for an 
energy window containg eleven eigenstates centered on the ll'th 
Bohr level. The thick dark line is the superposition of each en- 
ergy shell in the plane (periodic orbit) for energies in the range 
considered; the middle one corresponding to the Bohr-quantization 
energy. 

n = 304 n = 305 
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FIG. 4. Some of the individual Wigner functions on the (g2>P2) 
plane that were averaged in Fig.^. The plots are globally normal- 
ized to point out the enhanced amplitude of the scar for the n = 305 
state. The eigenenergy of this state is at distance ss 1.5A from the 
Bohr energy (A the mean level spacing in the chosen average win- 
dow). The black ellipse, in each graph, is the periodic orbit for the 
corresponding eigenenergy. 
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FIG. 5. The application of a particular nonlinear canonical 
transformation to the Hamiltonian (^) deforms the invariant plane 
qi = pi = into the surface displayed in the graph a, that also has 
pi = 0. The curve over this surface is the distortion of the original 
periodic orbit (ellipse displayed in the invariant plane). The graph 
b shows the new central surface for the deformed periodic orbit. 
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FIG. 6. Two of the many symmetric periodic orbits that have 
chords centered on the (92, P2) plane. The graphs are in position 
space, so that the momenta for the symmetric endpoints are dis- 
played as vectors. All centers have 52 ~ 0. The thin line is an 
cquipotential of (H) for the energy of the ll'th Bohr level. 
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